Ferromagnetic fluctuation and possible triplet superconductivity in Na^CoC^ • yH 2 0: 
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Spin and charge fluctuations and superconductivity in a recently discovered superconductor 
Na;cCo02 • J/H2O are studied based on a multi-orbital Hubbard model. Tight-binding parame- 
ters are determined to reproduce the LDA band dispersions with the Fermi surfaces, which consist 
of a large cylindrical one around the F-point and six hole pockets near the K-points. By applying the 
fluctuation-exchange (FLEX) approximation, we show that the Hund's-rule coupling between the Co 
t2 g orbitals causes ferromagnetic (FM) spin fluctuation. Triplet f v ( y 2_ 3x 2\-wsve andp-wave pairings 
are favored by this FM fluctuation on the hole-pocket band. We propose that, in Na^OoC^ • 1/H2O, 
the Co tig orbitals and inter-orbital Hund's-rule coupling play important roles on the triplet pairing, 
and this compound can be a first example of the triplet superconductor in which the orbital degrees 
of freedom play substantial roles. 
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PACS numbers: 74.20.Mn 

Recently, the first Co-based oxide superconductor 
Na^CoC^ • J/H2O with T c ~5 K was discovered by Takada 
etal. In this compound, edge shared CoOg octahedra 
form C0O2 planes similar to the high T c cuprates with 
Cu02 planes. However, the Co ions form a triangular 
lattice in the C0O2 plane. Mechanisms and nature of 
the superconductivity are currently attracting great in- 
terest because of a possible realization of unconventional 
pairing or resonating- valence bond (RVB) state 0, 
Although there are still some controversies, extensive Co 
NMR/NQR experiments have shown evidences for 11011- 
s-wave pairing 0, H Various experiments on the 
non-hydrate compounds show characteristic behaviors of 
strongly correlated electron systems Hp and in-plane 
ferromagnetic (FM) spin fluctuations 0, E3 • I n spite of 
these extensive studies, however, the pairing symmetry 
has not been determined. 

According to the triangular crystal symmetry, super- 
conducting states are classified as s-, p-, d-, f x (x 2 -3y 2 ) ( or 
next-nearest-neighbor /-), f y (y 2 -3 X 2 )j an d «-wave pair- 
ings. Previous microscopic theories using various model 
Hamiltonians have shown that various symmetries arise 
depending on the model. For example, if one uses an 
extended Hubbard model which has a charge-ordering 
instability in large- y region 0, 0], one obtains the 
next-nearest-neighbor /-wave pairing EH- Mean-field 
analyses of t-J model predict the chiral d xy + id 
wave pairing El El III El , and d xy ~t~ id x 2^y 
/-wave states are obtained in the single-band Hubbard 
model [h3 . I20I l2l| . We think that this variety of super- 
conducting states is due to frustration which does not 
give any specific momentum dependence of x«(?)- 

The above theories, however, assumed single-band 
models. As shown by LDA calculations |22tl23|. there is 
apparently an orbital degeneracy in Naa;Co02 systems. 
Actually two bands constructed from the three Co ti g or- 
bitals intersect the Fermi level. Furthermore, the electron 
fillings of these two bands are far from half filling. This 
means that multi bands (in other words, multi orbitals) 
contribute to low-energy electronic structure. This is in 



x 2 +y 2 ~ 

and 



marked contrast with the high-T c cuprates. 

In this letter, we study the electronic structure and su- 
perconductivity in Naa;Co02 • 2/H2O on the basis of the 
multi-orbital Hubbard model, which includes the Co ti g 
orbitals and the inter- and infra-orbital interactions. We 
find that several important and interesting aspects ap- 
pear which were not expected in the single-band model. 
One of them is a FM spin fluctuation which appears 
from the six hole pockets consisting of e' g orbitals as de- 
scribed below, and which is enhanced by the inter-orbital 
Hund's-rule coupling. Actually, an evidence for the FM 
fluctuation has been observed by Co-NQR experiments 
in several groups 0,H0l an d expected in the LSD A cal- 
culation [23] • We will show that this FM spin fluctuation 
leads to triplet pairing mainly on the hole pockets. Here, 
the disconnectivity of the Fermi surface plays an impor- 
tant role. First, we will deduce the tight-binding param- 
eters which reproduce the LDA band dispersions |2^ . l23| . 
Then we apply the fluctuation-exchange (FLEX) approx- 
imation to this multi-orbital Hubbard model. Generally 
speaking, the FLEX calculation is appropriate in the case 
near a quantum critical point because the critical en- 
hancement of fluctuation can be taken into account. We 
point out that Naa;Co02 ■ J/H2O can provide a first ex- 
ample of the triplet superconductivity which clarifies the 
important roles of orbitals on the superconductivity in 
the strongly correlated electron systems. 



We study the multi-orbital Hamiltonian given by H 
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kin. 
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and H int . = Hu + Hu>+Hj b +Hji, where d ima (d\ ma ) is 
the annihilation (creation) operator of an electron with 
spin cr(=1", J.) in the orbital m on the Co site i. Here, m 
runs over the xy, yz, and zx orbitals. The first term if cry . 
expresses the crystal field from the O ions acting on the 
Co tig orbitals. In Naj;Co02 ■ J/H2O, the CoOg octahe- 
dron is contracted along the c-axis. The ligand oxygens 
on the C0O6 octahedron with this distortion generate a 
trigonal crystal field, and the matrix element of D mn is 
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FIG. 1: (a) Triangular lattice and its Bravis vectors, (b) Band 
dispersions and (c) Fermi surfaces obtained from Hkin. +#cry., 
which reproduce very well the characters of the LDA re- 
sults |22j especially near the Fermi level. The bands far below 
the Fermi level are relatively different from LDA results, im- 
plying a possible importance of further long-range hoppings 
and the O 2p components. 

given by -j(5 m „ — 1). This crystal field lifts the local 
tig degeneracy into lower a\ g and higher e' g levels with 
energy splitting of A. 

The second term J/kin. represents the kinetic-energy. 
Here, \i is the chemical potential, and t^ n is the hop- 
ping integral between the m orbital on the i site and 
the n orbital on the j site. In order to reproduce the 
LDA band structures, we find that up to the third 
nearest-neighbor hoppings between orbitals are neces- 
sary. In the /c-space, flkin. is rewritten as J/kin. = 

Ek,m,„, CT £ r4m^kna - fi £k,m,<7 4m<Ama with 

e 77 = 2ti cos kZ< + 2i 2 [cos fc 77 + cos^ 7 + fc 77 )] (1) 

+ 2i 4 [cos(2fc 77 + fc 77 ) + cos(fc 77 - fc 77 )] 

+ 2i 5 cos2fc; 77 , 

e^' = 2i 3 cosfc 77 ' + 2i 6 cos2fc 77 ' (2) 

+ 2t 7 cos(^ 7 ' + 2fc 77 ') 

+ 2i 8 cos(fc 77 ' - fc 77 ') + 2i 9 cos(2fc 77 ' + fc 77 '). 
Here, 7 and 7' represent xy, yz and zx orbitals and 

k xy,xy = k xy,zx = ^ fe »»,x» = fc x»,.» = ^ k y Z ,y Z = 

h*™* = fc 2 , kl Z ' VZ = k* y ' VZ = -(kl + h 2 ), = 
k yz,zx = _( ki + fe2 ) and f.zx,zx = k yz,zx = ^ respeo _ 

tively. ki and fc 2 are the components of the wave vector 
k of the triangular lattice spanned by d[ and 02, respec- 
tively (see Fig. [T] (a)). In Fig. ^ wc show the band dis- 
persions (b) and the Fermi surfaces (c) obtained from 
#kin. + #cry. in the case of t\ = 0.45, ti = 0.05, t 3 = 1, 
U = 0.2, t 5 = -0.15, t 6 = -0.05, t 7 = 0.12, i 8 = 0.12, 
< 9 = -0.45 and A = 0.4. Both reproduce well the LDA 
results 2 2j , particularly near the Fermi level. In the fol- 
lowing, we use these parameters and i 3 = 1 as the energy 
unit. 

As shown in Fig.Q] the Fermi surfaces consist of a large 
cylindrical one around the T-point and six hole pockets 
near the K-points. The cylindrical one has a dominant 



ai g -orbital character, while the six hole pockets have an 
e'g-orbital character. Koshibae and Maekawa j2|| dis- 
cussed hidden Kagome lattices using the second-nearest 
neighbor hoppings. However, in their case, the obtained 
Fermi surfaces have six large hole pockets around the 
M-points in contrast to the LDA result. In this com- 
pound, the hybridizations between the neighboring Co 
tig and e g orbitals are relatively strong because of the 
lattice structure with edge-shared CoC-6 octahedra, and 
a hopping process via the unoccupied e g orbitals aids the 
third nearest neighbor hopping but not the second near- 
est one. Thus, it is essential to take into account up to 
the third nearest-neighbor hoppings for reproducing the 
LDA result. In Fig.U(b), the bands far below the Fermi 
level are relatively different from LDA results, implying 
a possible importance of further long-range hoppings and 
the O 2p components. However, the following results do 
not depend on these details far below the Fermi level. 

The last term in H, H lnt , represents the on-site d-d 
Coulomb interactions, where Hjj and Hu> are the intra- 
and inter-orbital Coulomb interactions, respectively, and 
iJ/ H and Hji are the Hund's-rule coupling and the pair- 
hopping interactions, respectively. These interactions are 
expressed using Kanamori parameters, U, U', Jh and J', 
which satisfy the relations; U' = U — 2Jh and Jh = 
J'. The value of U has been estimated as 3-5.5 eV by 
the photoemission spectroscopy and as ~3.7 eV by 
the ab-initio calculation [2^]. The value of Jh for the 
Co 3+ ion is 0.84 eV. Thus, the ratio Jh_/U, which gives 
the strength of Hund's-rule coupling, is 0.15-0.28 in this 
compound. 

We analyze this multi-orbital Hubbard model by ex- 
tending the FLEX approximation to the multi-orbital 
case (23. In the FLEX approximation, RPA-type bubble 
diagrams and ladder diagrams are taken into account. By 
determining self-consistently the renormalized Green's 
function and fluctuation exchange self-energy, the effects 
of mode-mode coupling between charge and spin fluctu- 
ations and quasi-particle life time due to damping are 
incorporated, which weaken the Stoner-type instability 
expected in the mean-field type analysis. In the present 
three-orbital case, the irreducible Green's function G, the 
non-interacting Green's function G^ and the self-energy 
E are expressed in the 3 x 3-matrix form corresponding to 
the xy, yz and zx orbitals. The irreducible susceptibil- 
ity x° and the effective interaction V have a 9x9-matrix 
form. The Green's function satisfies the Dyson-Gor'kov 

equation G(/c) _1 = G ( ^(k) -1 - ±(k) with 

S TOn (fe) = — y] y] V^ m ,un(Q)Gfiu{k - q)- (3) 

q fiV 

Here, the matrix element of the effective interaction is 
given by 

V m n,nv(q) = [^U s X s (q)U s + \u c X c (q)U c 

-i(L> s + U c ) X °(q)(U s + u c ) + \u s - \u c ] mn ,,M) 
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K r M K largest A: (a) /j / (j / 2_ 3a .2)-wave gap and (b) p y -wave gap. 



FIG. 2: (a) Maximum eigenvalue of x s showing the evolution 
of FM fluctuation, as the strength of Hund's-rule coupling 
Jh/U increases, (b) Fermi surfaces for interacting case with 
[7=6.0 and J H /(7=0.18. 

with x s (q) = [I - X {d)U s ]- l x\d) and x c (q) = [I + 
X°(9)?7 C ] _1 X (9)- The irreducible susceptibility is given 
by 

X%m,v,M) = -Jf G A»»( fc + q)G m {k). (5) 

k 

In the above, T is temperature, k = (k, iuo n ) and q = 
(q, ivi) with uj n = (2n — l)nT and v\ — 2lnT, N is the 
number of sites, and J is a unit matrix. The interaction 
matrices U s and U c are represented as 



u s = 


"t/f 
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u c 2 . 



(6) 



where [Ul] aa ,bb {[U{} aa ,bb) is U (U) for a = b, and J H 
(2C/' — Jh) for a ^ b. On the other hand, [£/|] a &,cc2 
([C/llab.cd) with a ^ b and c ^ rf is given by /7' 
(-[/' + 2 J H ) for (a, b) = (c, d) and is given by J' (J') 
for (a, 6) = (rf, c), and is for the other cases. 

To discuss the nature of superconductivity, we solve 
the following Eliashberg equation; 

Wmn(k) = ~EEE r l^(9)^(^-?) 

q a/3 liv 

x G atl {k-q)Gf 3u {q-k). (7) 

The singlet and triplet pairing interactions T s and T* 
have 9x9- matrix forms as, 

f s (9) = \u s x s {q)u a - \u c x c {q)u c + \{u s + u c ), (8) 

f*(g) = - l -U s x s {q)U s - \u c X c {q)U c + \{U S + U% (9) 

The eigenvalue A is a measure for the dominant super- 
conducting symmetry, and becomes unity at T = T c . 

Calculations are numerically carried out with 32x32 
/c-meshes in the first Brillouin zone, and 4096 Matsubara 
frequencies. The value of U is fixed at 6.0, and tempera- 
ture T is fixed at 0.02. First, let us discuss the spin and 
charge fluctuations. Figure [21 (a) shows the maximum 
eigenvalue of x s obtained by solving the above FLEX 



equations self-consistently for various values of Jh/U at 
T =0.02. Apparently, a peak structure around T-point 
is strongly enhanced as Jn/U increases, while the small 
peak at the M-point hardly changes. Here, we also show 
the Fermi surfaces for interacting case with U =6.0 and 
Jn/U=0.18 obtained by the present FLEX calculation 
in Fig. [3(b). The peak at the M-point originates from 
the inter-hole-pocket scatterings with wave numbers of 
qi and q 2 shown in this figure |3l|. Since the value 
of Jh/U is 0.15-0.28 in the actual compound, the en- 
hanced peak at T-point at Jb_/U =0.18 indicates that 
the FM fluctuation is dominant in this compound. As 
mentioned before, FM fluctuation has been observed ex- 
perimentally [j, |24| . On the other hand, x° does not 
exhibit any remarkable structures, which implies an ab- 
sence of charge- and orbital-density-wave instabilities. 

By solving the Eliashberg equations, we find that 
triplet /j / (j / 2_ 3a ,2- ) -wave or p-wave pairings have the largest 
A in the calculated parameter region of < Jh/U < 0.30. 
Although the values of A are nearly the same, A for 
/j / (y2_ 32 ,2)-wave state is slightly larger in the region of 
Jh/U < 0.2 while A for p-wave state is larger in Jn/U > 

0. 2. The other pairing instabilities are considerably weak. 
We have also studied a second-order perturbation theory 
and have obtained consistent results [3fj|. 

Figure[3shows the fc-dependence of obtained supercon- 
ducting gaps with (a) / y ( y 2_ 3:z .2- ) -wave and (b) p y -wave 
symmetries. As for the p-wave, the p x and p y states 
are degenerate in the triangular lattice. In addition to it, 
there are three choices of the direction of d- vector (x,y,z), 

1. e., the S'z-component of 15=1 Cooper pair. Below T c , a 
linear combination of these six states should be realized, 
namely, i(p x ±ip y ), xp x ±yp y and ip y ±yp x . All of their 
gaps are represented as yj A x (k) 2 + A y (k) 2 with A x (k) 
(A y (k)) being an order parameter for p x (p y ) state. This 
gap structure has line nodes on the a\ g Fermi surface 
since the line nodes of p x and p y states have nearly six- 
fold symmetry as shown in Fig. 01(b). It is apparent that 
the magnitude of gap is large on the e' g band with six 
hole-pocket Fermi surfaces, while it is considerably small 
on the dig band. The geometry of the Fermi surfaces 
plays a very important role: Because of the disconnec- 
tivity of the Fermi surfaces, the nodal lines run between 
the e'g hole pockets and do not intersect them. Conse- 
quently, the gap is fully opened on each e' g Fermi sur- 
face, while a\g Fermi surface has line nodes. As discussed 
before the discovery of Na^Co02 • yH^O, this kind of 
disconnectivity favors superconductivity by avoiding the 
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disadvantage from nodes on the Fermi surface |32(. The 
dominant contribution of the e' g band to the supercon- 
ductivity is also attributed to the van Hove singularity 
(vHS) in the e' g band. As can be seen in the band struc- 
ture in Fig.^ there exist saddle points near the K-points. 
This leads to a large density of states (DOS) near the 
Fermi level, which was actually reproduced in the LDA 
calculation p^ . 

Experimental estimate of the value of Jn/U has an 
ambiguity and is ranged from 0.15 to 0.28. Thus, both 
fy(y 2 -3x 2 ) state (Jft/U < 0.2) and p state {J-r/U > 0.2) 
can be candidates for the pairing state in the actual ma- 
terial. Now, let us discuss the present results with ex- 
periments. Both /-wave and p-wave gap structures ob- 
tained here are consistent with several NQR/NMR ex- 
periments 0, |(J, which suggest existence of line nodes. 
Time-reversal symmetry observed in //SR [3j| excludes 
the possibility of chiral p-wave state i{p x Hp y ). Further- 
more, the recent measurements of Knight shift |2J, y4j 
and H c2 [l^, |3(| are consistent if the d- vector is fixed to 
be parallel to the plane. This supports the xp x ± yp y 
or xp y ± yp x states. However, a contradictory result 
about the constant Knight shift below T c has also been 
reported 3 so that further careful study is needed. In 
principle, the spin-orbit interaction lifts the degeneracy 
and determines which state is stabilized. This is an in- 
teresting future problem. 

The existence of e' g hole pockets is still controver- 
sial. Several band calculations have predicted it [22L |23| 
while a recent LSDA+U calculation showed only a large 
di s Fermi surface [13 . However, it was found that the 
L(S)DA+U cannot reproduce the experimental optical 
spectra [2^ , suggesting that this method is probably ill- 



suited for the present Co oxides. The angle-resolved pho- 
toemission spectroscopy has been performed only for the 
non-hydrate material || and the hole pockets have not 
been detected, which may be due to a surface effect or 
some matrix elements. Our calculation shows that the 
hole pockets play a substantial role for superconductivity. 
Indeed, if we study the situation without hole pockets by 
tuning the transfer parameters, we find that the pairing 
instability is strongly suppressed and the realization of 
superconductivity with T c ~ 5 K is difficult. Our results 
strongly suggest the existence of hole pockets. 

In summary, we have studied the spin and charge fluc- 
tuations and the superconductivity in Naa;Co02T.3H20 
on the basis of the FLEX analysis of the multi-orbital 
Hubbard model. The Hund's-rule coupling between 
the Co t2 g orbitals gives rise to the FM spin fluctua- 
tion in agreement with the NQR results from several 
groups |1 H ||. The triplet pairing instability medi- 
ated by this FM fluctuation arises mainly on the e' g band 
with the pocket Fermi surfaces. The hole-pocket Fermi 
surfaces as well as a large DOS due to the vHS near the 
Fermi level are crucially important for the superconduc- 
tivity. We have pointed out that this material can be 
a first example of the inter-orbital-interaction-induced 
triplet superconductivity. This mechanism is a univer- 
sal one which can be expected in superconductors with 
multi-orbital degrees of freedom. 
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